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! Abstract 

' A system of electrons in the two-dimensional honeycomb lattice 

with Coulomb interactions is described by a renormalizable quantum 
field theory similar but not equal to QED3. Renormalization group 
techniques are used to investigate the infrared behavior of the system 
that flows to a fixed point with non-Fermi liquid characteristics. There 
are anomalous dimensions in the fermionic observables, no quasiparti- 
^ ' cle pole, and anomalous screening of the Coulomb interaction. These 

' results are robust as the Fermi level is not changed by the interaction. 

The system resembles in the infrared the one-dimensional Luttinger 
liquid. 
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1 Introduction 



This paper addresses the problem of interacting electrons in a band of a two- 
dimensional crystal. This problem uncovers all the macroscopic properties 
of the materials of which the maximum interest is centered at the present 
on the high- Tc superconductors, in particular, on the nature of their normal 
state. To fix ideas and some notation let us review the highlights of the 
problem. Start with a free theory of, say, electrons in a metal, i. e., a free 
Fermi gas whose ground state is known and whose low energy elementary 
excitations are electrons and holes obeying a free dispersion relation. Then 
switch on an interaction among the electrons and ask the question of what 
the new ground state is and what is the nature of the low energy excitations 
of the full interacting theory. This program has been successfully addressed 
by Landau in three spatial dimensions d where it leads to the famous theory 
of Landau's Fermi liquid [IH and the concept of quasiparticles. Things are 
different m d = 1 where an exact solution can be found and the behavior of 
the system is always of the Luttinger typepf. 

The case c? = 2 is physically very interesting both in practice (most of 
the high-Tc materials are organized in two dimensional layers and experi- 
mental data are available) and under a theoretical point of view. Despite the 
great effort devoted to it, no exact solution is known at the present and the 
situation remains controversial]^. The ultimate origin of the trouble lies on 
the very complex nature of the many-body vacuum and the corresponding 
difficulties of doing perturbations around the Fermi surface. 

The theory of the Fermi liquid and the perturbations of it that lead to 
the BCS and charge density wave instabilities have recently been rephrased 
as effective field theories 0, ||, ^ to which renormahzation group techniques 
can be applied 0]. The main idea of this approach is that, although we are 
dealing with problems that naturally have an ultraviolet cutoff po (the width 
of the band), we are only interested in the physics at a scale much lower 
than it, namely, a very narrow strip above and below the Fermi surface. By 
integrating out in the path integral all the momenta above a certain cutoff 
P < Po, we end up with an effective field theory whose couplings are cutoff 
dependent. The renormahzation group transformations are rescalings of p. 
In this way one can visualize the Fermi liquid as a fixed point of this RG 
transformations (it is gapless) that fiows towards a BCS state when a weak 
attraction is introduced between the quasiparticles. 

When comparing the effective field theories that arise in this context with 

■^Onc of the clearest examples of the mentioned difhculties is the "infrared catastrophe" 
discovered by Anderson that is attracting a major interest in relation with the study 
of the high Tc superconductivity. 
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standard quantum field theory, we encounter two main differences. The first 
one is the nonrelativistic nature of the many body physics where there is 
a parameter, the Fermi velocity, whose typical value is much smaller than 
the speed of light. This prevents the use of the general theorems concerning 
renormalizability. The more severe difference lies in the above mentioned 
nature of the many body vacuum. When treating the interacting system as 
a perturbation of the "free" system constituted by the fixed point, pertur- 
bations are to be described around a Fermi surface in ci = 3, a Fermi line in 
= 2, and Fermi points in d = 1. We are then faced with couplings that 
are functions of the shape of the Fermi surface which is also modified by the 
interaction. Only in the d = 1 case we deal, as in quantum field theory, with 
a true vacuum around which one can expand and get a bunch of coupling 
constants. 

In ref. we developed a model for the electronic structure of free 
electrons in a two-dimensional honeycomb lattice. This lattice constitutes the 
basic building block of periodic systems made of identical atoms with planar, 
threefold coordination. That is the case of graphite, built of carbon atoms. 
The rapidly growing family of fullerene compounds |^ |10| can be thought 
of as closed, curved surfaces derived from the planar honeycomb lattice. 
These compounds show a variety of unexpected features, the most striking 



of them being the superconductivity of doped crystals of Cgo molecules [|TT 
Moreover, this superconductivity can be considered a form of high-T^ 
superconductivity, as the number of electrons available to form Cooper pairs 
is unusually low. 

The main feature that makes the honeycomb lattice the object of our 
attention in the present context is that its Fermi level at half-filling consists of 
exactly two Fermi points. This peculiarity opens the possibility of performing 
a complete analysis of the problems mentioned above in analogy with the 
d = 1 case. In particular it makes of this lattice an ideal laboratory for 
investigating departures from Fermi liquid behavior in c? = 2 which have 
shown very ellusive in more conventional approaches. The major interest of 
this investigation relies on the fact, made clear by now from the experiments, 
that the normal state of high-Tc materials is not of Fermi liquid type. Besides, 
it can be directly applied to the physics of the fuUerenes that, as mentioned 
above, are two-dimensional graphitic structures at half-filling |TE|. 

The paper is organized as follows. In section 2 we review results on the free 
electrons on the honeycomb lattice that lead to the low energy effective action 
used as a starting point of our analysis. We emphasize the point that, unlike 
what usually occurs in many body physics, our effective action constitutes a 
genuine quantum field theory with a single, well-defined, vacuum state. Next 
the free theory is completed with the addition of the Coulomb interaction. In 
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the field theory, this is mediated by an electromagnetic field whose propagator 
reproduces, in the nonrelativistic limit, the usual four-fermion interaction 
considered in the many-body approach. We end up with a standard local 
quantum field theory. A discussion is done on the differences to be expected 
by the use of the relativistic versus the instantaneous electromagnetic field 
propagator illustrated by the computation of the density of states in both 
cases. The results announce an infrared instability which is clarified in section 
4 with the help of the renormalization group. 

Section 3 is devoted to the study of the renormalizability of the model. It 
is shown that, despite its nonrelativistic nature and due to its gauge invari- 
ance, it is renormalizable. The Fermi velocity and the electron wave function 
get renormalized while the electric charge does not. 

Sections 4 and 5 contain the analysis of the RG flow of the model and 
its physical implications. The first section is devoted to the nonrelativistic 
limit of the theory. We see that the perturbative expansion has an effective 
coupling given by the ratio of the square of the electric charge and the Fermi 
velocity. The RG shows that this coupling decreases to zero in the infrared 
although there is no fixed point in the nonrelativistic regime. We then in- 
vestigate the possibility of having a natural length scale in the theory and 
make some comments on the system away from half-filling. Finally we find 
an anomalous dimension for the electron wave function at two-loop order 
that points towards a non-Fermi liquid behavior of the system. 

The previous analysis is completed in section 5 where we study the rela- 
tivistic regime which is the relevant regime in the infrared. We find a fixed 
point of the RG fiow obtained when the Fermi velocity equals the speed of 
light. Next we compare our model with standard RG analysis of the Fermi 
liquid theory and find significant differences both from the technical and 
from the physical point of view. The most significant one is the anomalous 
dimension found for the electron Green function. We discuss the non renor- 
malization of the electric charge and establish the non-Fermi liquid nature 
of our fixed point. 

We set our main conclusions and a summary of the more relevant points 
of the article in section 6. 

2 Many-body theory of the 2D layer 

In this section we first review the one-particle electronic states in the two- 
dimensional honeycomb lattice, which serve afterwards as a starting point 
for the many-body problem. The present one-particle description has been 
worked out with complete detail in ref. [§, regarding the conduction band of 
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a two-dimensional graphite layer. The honeycomb lattice has two different 
atoms per primitive cell, so that two degrees of freedom arise in a variational 
computation of the energy eigenstates, corresponding to the respective or- 
bitals for the two atoms. In the tight-binding approximation the problem 
reduces to the diagonalization of the one-particle hamiltonian 

H=-tY: at a, (1) 

<i,j> 

where the sum is over pairs of nearest neighbors atoms i,j on the lattice and 
di , dj are canonically anticommuting operators 

{ai, aj} = {a+, a+} = {a^, a+} = 6ij (2) 

The states which are simultaneously eigenvectors of the hamiltonian and of 
the lattice generators have to be of the form 

^ = E c.e'^'"'4 \o) + E coe'^'^'at \o) (3) 

i • i o 

with coefficients c, and Co for black and blank points, respectively, as depicted 
in figure 1. Obviously, under the action of black points are mapped into 
blank points and viceversa. We have, indeed, 

= -^E E c.e'^"'at\0) - tE E c.e'^''^at\0) 

i • <i,j> i o <i,j> 



where {uj} is a triad of vectors connecting a • atom with its nearest neigh- 
bors, and {vj} the triad made of their respective opposites (see figure 1). 
Therefore, the state is an eigenvector of H provided that the coefficients 
c, and Cq cirG solutions of the eigenvalue problem 



-tEje'^^"^^ \ ( c, \ ( c, 

-,zk-v 



^(k) • (5) 



The parameter t measures the hopping between nearest neighbors and pro- 
vides the natural energy scale in the one-particle description. The diagonal- 
ization of gives the band of levels 



E(k) = ±ty 1 -|- Acos'^—^akx + Acos—^ak^ cos-aky (6) 
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the parameter a being the lattice spacing. We will be interested in the 
consideration of this band at half-filling, which is the pertinent instance for 
the carbon-based materials and, in particular, graphite layers and undoped 
fullerenes. The dispersion relation has the amazing property that the 
Fermi level at -E'(k) = is only reached by six isolated points in momentum 
space (see figure 2). These are the six corners of the first Brillouin zone 
and, due to the periodicity of the dual lattice, they correspond to only two 
independent states. The existence of a finite number of Fermi points is quite 
unusual in more than one spatial dimension. It has important consequences 
since it makes possible to encode the low-lying excitations about the Fermi 
level into a simple field theory. When considering an arbitrarily large sample, 
it is appropriate to amplify the region about any of the two independent Fermi 
points by taking the limit a — 0. In either case we obtain for the operator 
H in equation (|) 

n = ha (T-6k + O (^{a Sk f) (7) 

where dk is the displacement about the Fermi point and cri,a2 stand for 
the two Pauli matrices. The conclusion is that, in the long- wavelength limit 
a — s> 0, the one-particle (energy) eigenstates of the electrons in the honeycomb 
lattice are given by the field theory of two massless Dirac spinors (one for 
each independent Fermi point) in two spatial dimensions^. This is all about 
the free theory. 

We now turn to introduce the electronic interaction by applying the 
method of second quantization. We remark that, proceeding in this way 
with the above set of one-particle states, we end up with a formalism which 
falls more into the framework of quantum field theory than into that of quan- 
tum statistical theory. The difference is subtle but significant. In quantum 
statistical theory the ground state of the system is, in general, a condensate 
with an extended Fermi surface (in more than one spatial dimension) and 
the elementary excitations (quasiparticles) are only stable right at the Fermi 
level. The second quantized theory built from the dispersion relation (|^) 
has, however, a definite vacuum in the sense of quantum field theory. Its 
excitations are similar to the particle-hole excitations of the Dirac sea. We 
may consider separately, for instance, each of the two Dirac spinors and work 
out the many-body electron propagator of the free theory 

GW(t,r;t',r') = -z{T^ i{t,r)^j{t' ,r')) (8) 

We denote by the electron operator \1/ in the interaction representation 
while = \1/|(T3. In terms of the respective modes V'k,+ , ^k,- for unoccupied 

"^Similar dispersion relations can be found in higher dimensions. The best known of 
them is the one derived from the 3D diamond structure with one state per site|L4|. 
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and occupied levels we liaveQ 

^i{t,r) = }_^e e ' ' V'k,+ak,+ + 2^e e ' ' V'k-ak- (9) 

k k 

where v = 3ta/2 is the Fermi velocity and the operators ak,± satisfy anti- 
commutation relations 

{ak,±,a+_i} = 5(k-k') (10) 

Taking the statistical average with the system at half-filling, we get for the 
Fourier transform of the propagator 

G'°>(..fc)^. -'r+y'' (11) 

{70,7} is a standard set of 7-matrices satisfying 

{70,70} = -2 {7.,7j} = 25,,, {70,7.} = (12) 

and related to the above quoted Pauli matrices by (73 = —i-yo, o\ = icr^-yi, 02 = 
zcr372. The above expression coincides precisely with the Feynman propaga- 
tor for a massless spinor in 2 + 1 dimensions, which confirms our assertion 
that the system at half-filling is equivalent to the Dirac sea in quantum field 
theory. We must notice however the presence of the Fermi velocity v in front 
of the spatial part of the scalar product, which is a feature of the many-body 
theory. 

In the nonrelativistic approximation, the introduction of the electronic 
interaction leads, in first instance, to the second quantized hamiltonian 



S r 

Hcouiomb = -tay rfV ^(r)7 • V^(r) 



+ [£n /,v/(^i)^3^(rO^(r2V3^(r2) 3 
2 J J 47r |ri — 

We attach to this hamiltonian the name of "Coulomb" since it describes the 
instantaneous interaction between electric charges. The model given in terms 
of the Coulomb interaction is, however, a highly nonlocal field theory, what 
makes very awkward its investigation from a formal point of view. Both 
for the sake of studying the properties of the quantum theory as well as for 
computational purposes, a complete description of the interaction with the 
electromagnetic field is more desirable. One of the points that we support 

^From now on we replace 5k by the wavevector k with origin at the Fermi point. 
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in what follows is that, actually, the model of interacting electrons is sensitive 
to retardation effects of the electromagnetic propagation. The reason for such 
unconventional behaviour is the massless character of the spinor This 
property is not accidental in the description of the two-dimensional layer, 
as long as the expansion (|^ gives only one marginal operator dictating the 
long-distance behaviour of the free theory. 

We propose, therefore, a quantum field theory description based on the 
second quantized hamiltonian[| 

H = ^taj d\ ^(r)7 • V^(r) - e j d\ j^A^" (14) 

The interaction of the electromagnetic field and the electrons in the layer 
is described in the standard fashion, by coupling to the conserved current 

~ (i 1^70^, i t;'^7^) (15) 

This poses some technical problems since the electromagnetic field propagates 
in three-dimensional space while we want the dynamics of the electrons to be 
confined to the two-dimensional layer. Although this may not be achieved in 



general, it turns out to be possible by specializing to the Feynman gauge [|T5 
which enforces the constraint 

V^A^ = (16) 

and has also the property of placing the A^ field in the same direction as the 
electronic current. In this gauge we have the propagator 

,4, „zk-(r— r') p—iuj{t—t') 
{TA.it, r)A^it', r')) = -^6,. [ -— . (17) 

and the coupling to the 2 + 1 dimensional current is perfectly consistent. In 
the computation of electronic properties we have just to take care of placing 
the points r, r' on the plane of the two-dimensional layer. With this descrip- 
tion we certainly get more than we would need in a nonrelativistic theory, 
but we can recover at any time the nonrelativistic limit of all quantities by 
expanding in powers of v/c. 

It is important to stress that we have to expect, in general, different an- 
swers in the computation of a given observable taking (|T^) or (|T^ as starting 
point. The dimensionality of the system is low enough so that the massless 
condition of the spinor reflects in the appearance of infrared divergences in 
perturbation theory. Thus, although we should not expect any discrepancy 

^Unless otherwise stated, we work henceforth in units Ti = c = 1. 
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in the computation of local quantities (like cutoff dependent quantities) in 
the nonrelativistic model (p!3D and in the model given by (|1^ in the limit 



f/c — s> 0, the response of nonlocal quantities to the interaction turns out to 
be quite different in the infrared regime. From a technical point of view, per- 
turbation theory looses predictive power because of the mentioned infrared 
instabilities, and one has to resort to more sophisticated methods in order to 
obtain information from the quantum theory. 

Let us illustrate the above fact with the computation of the quantum 
corrections to the density of states n{uj) near the Fermi level {u = 0). This 
observable is given by 



n(LU\ 



Im J d^k Tr [G(cj, k)a^] (18) 



In the theory with the instantaneous Coulomb interaction, the density of 
states is not modified to first order in perturbation theory. The instantaneous 
interaction can be obtained from (|1^, for instance, by considering formally 
its expression for c oo. In this limit we have, computing always at points 
r, r' in the two-dimensional layer 



{TA^{t,v)A,{t\v')) ^ -z6,,J 



d'^kduj rdk.e"^'^' ' ' e 



iik-(r-r') -iuj{t-t') 



(27r)3 J 27r + kl 



le 



Then, it becomes obvious that this kind of interaction cannot modify the 
uj dependence of the propagator G, to the first perturbative order (see also 
the discussion at the end of section 4). However, as mentioned above, the 



propagator (|T^) cannot be safely used in the computation of nonlocal quan- 
tities. This remark is appropriate to the case of the density of states, whose 
determination in the quantum field theory given by (|^) goes as follows. 

The inverse propagator may be decomposed, as usual, into a free part 
and the electron self-energy S(c(j, fc) 

5 = ^-2 (20) 

S(ct;, k) has a perturbative expansion in terms of one-particle-irreducible 
diagrams [|l|, which are in general ultraviolet divergent. Since our model can 
be treated as a genuine quantum field theory we expect these divergences 
to be local and susceptible of being absorbed into 1/G^^\ For the moment 
we are only interested in finite corrections and will leave the issue of renor- 
malization to the next section. To first order in perturbation theory, the 
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relevant diagram is shown in figure 3 and gives (in the hmit v/c 0) the 
finite contribution to S 

Eji(uj,k) = i———v J • k / dx— 7-—:rlog\v'^k'^x-uj'^x(l-x)\+0(e'^) 

loTT^ JO [1 — X + v^xy ^ ' 

(21) 

When inserted in equation ([T8|) the above expression gives the first quantum 
corrections to the density of levels 

/ \ 2^1 1 , 2 12/"^, Vl - X 

n(w) = 27r — 1 + — /o^ ^ + — e / rfx- , ^ a2 ^ 

"^IGtt^ ^ (1 — X + f ^x)^ X ^ 

It is already appreciated that the perturbative expansion of n{uj) is not well- 
defined, since the last integral in (|22| ) is clearly divergent. Even if we could 
manage to regularize in some way this infrared divergence, we still would 
be left with a more serious problem. At the point near which we want to 
measure the density of levels, i.e. near a; = 0, the first perturbative correction 
dominates over the zeroth order term, and one can presume that higher order 
terms can produce increasingly large powers of log uj. Perturbation theory is 
clearly not reliable in the infrared regime, but the kind of infrared instability 
described points at new physics near the Fermi level. 



3 Renormalization 

Henceforth we will pay attention to the local quantum field theory (0), 
which offers the possibility of analyzing the infrared behaviour by applica- 
tion of renormalization group methods. We study first in this section the 
renormalization properties of the model. Though the many-body theory of 
the 2D layer has a natural ultraviolet cutoff, we insist in absorbing ultra- 
violet divergent contributions into bare parameters of the theory since this 
is a way of extracting relevant information about the scaling properties in 
the infrared|l^, As we have mentioned before, the renormalizability of a 
quantum many-body theory is, in general, a nontrivial issue, since the cou- 
plings in the effective theory may depend on the point chosen on the Fermi 
surface. In our case things become much simpler as our low-energy effective 
theory is a genuine quantum field theory. Regarding the usual considerations 
for renormalizability in field theory, the only condition that we are missing 
is relativistic invariance since the Fermi velocity of the electrons v = 3ta/2 
does not match the speed of light c that appears in the dispersion relation of 



9 



photons. For this reason it is still necessary to check that the ultraviolet di- 
vergences of the theory can be absorbed into a redefinition of the scale of the 
fields and of the parameters in the hamiltonian. We present here the analysis 
at the one-loop level, and conclude with a compelling argument showing that 
the model is renormalizable to all orders in perturbation theory. 

The conditions which make the theory renormalizable can be stated more 
clearly in the lagrangian formalism. Corresponding to the hamiltonian 
we have the action 

S = J dtd^r *'(-7o<9o + v-f • V)^ 

-ie J dtd^r "^{--foAo + v-f ■ A)^ (23) 

In the quantum theory the coefficients of the terms in the action have to be 
adjusted so as to cancel the ultraviolet divergences. We should start therefore 
with the bare action 

Share = Zkin j dtd^T ^(-70^0 + Z^V^i • V)^ 

-Z,nt ie j dtdh ■*(-7oAo + Z,v-i ■ A)^ (24) 

The renormalization coefficient Zkin may be determined from the electron 
self-energy E(a;,fc), as well as the renormalization coefficient for the Fermi 
velocity Z^. The coefficient Zint is determined from the renormalization of 
the time component of the interaction vertex. It turns out therefore that 
the spatial components have necessarily to be renormalized by the product 
ZintZv This is a nontrivial check of the renormalizability of the theory, which 
arises in the absence of relativistic invariance. Provided that it is fulfilled, the 
renormalized theory may be made finite in terms of renormalized quantities 

Vbare = Z^Vji 
ZeCR 



bare 



^bare = 

^bare = Zj A'^ (25) 

In the previous section we already quoted the one- loop renormalized con- 
tribution to the electron self-energy T,{uj,k). We deal here with the ultra- 
violet divergent part of this object, regularizing it by working in analytical 
continuation to dimension d = 3 — e[|r3|. From the diagram in figure 3 we 
have 
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= -ZvE,i(-7ou; + Z^v ^ • k) - i— -fou (1 - 2v^) / c/x- ^ — - — 

077^ ^ ^ JO 1 — X + f ^X £ 

— i — r V J • k / c/xt ^— TTT- + finite terms + O(e''026) 

Stt^ Jo (1 — X + v^x)^ e 

In order to ensure the finiteness of the electron propagator we take the renor- 
mahzation coefficients 

Z^ = 1 + (l - 2v') f dx 1 + O(e^) (27) 

svr"^ ^ ' JQ 1 — X + f ^x e 

Z„ = 1 / dx 



n"^ Jo (1 — X + v'^x)'^ e 

e^ / ^ 2\ /"S 1 



(l - 2v^) \ dx + O(e^) (28) 

87r"^ ^ ^ JO 1 — X + f ^x e 

These coefficients have a more comphcated structure than those of a rela- 
tivistic theory, as long as v does not match the speed of light — which, in 
our units, means that v ^ 1. Z^ has, for instance, an infinite power series 
expansion in v 

167r2 I w V2'2'2' / V J V2'2'2' 

+4 (l - 2v') F (l, 1; I v'^ - 4F (l, 2; | } ^ + O (e^) (29) 

We will later use this expression in the consideration of the nonrelativistic 
approximation v 0. 

Now we come to the renormalization of the interaction term in the action 
(^4|) . The ffist quantum corrections to the vertex are given by the diagram in 
figure 4, and we denote them symbolically by F^. The divergent contribution 
to the time component of the vertex is actually different than that to the 
spatial components. We find that 

/_ „ 2^ /"^ ' 1 — X 



lovr^ ^ ^ JO a/x (f ^ + (1 — v^)x) 

—z 7o [1 — 2v ) dx 



^vr2 ^ Jo ^ {v^ + {1 - v^)xf e 

+ finite terms + O(e^) (30) 

r = -i-r^v-f dx r— - + finite terms + 0{^}) 

lOTT^ JO a/X (f^ + (1 — f^jXj £ 

We may think of as a correction to the bare vertex in (EW)- We have then 

g2 \ /'-'^ 1 X 1 

^mt = 1 - (l - 2t;^) / dx 

lOTT^ ^ ^ JO a/x (i;^ + (1 — f jxj £ 
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By taking the quotient between (|33D and (|32D, we get an independent deter- 
mination of the renormahzation coefficient Zy 

Zy = 1 — -^—^ f dx- ^ 



IQtt'^ Jo y/x {v"^ + {1 — v^)x) \ v'^ + {l — v^)xj £ ~^ ^ 

= 1 _ ^ - 2v^)F (-, --; vA + -v' f 1 - 2v^) F (2,2; -; 

167r2 1 i;^ ^ V2' 2' 2' / 3 V ^ V 2' 

+2nvF (i I i; - 8v^F (l, 2; | } J + O (e^) (34) 

It is not obvious that this expression coincides with the alternative form (p9|), 
though it has to for the theory to be renormahzable. In fact, after using the 
recursion relations which hold for hypergeo metric functions [0] one can show 
that the two functions in (|3^ ) and (^) are identical. We then have the nice 
situation of dealing with a nonrelativistic quantum field theory in which the 
Fermi velocity of the electrons is renormalized exactly in the same way in 
the different terms of the action. 

This precise coincidence is not casual. One is led to suspect that it has to 
be consequence of some symmetry operating in the theory, which turns out 
to be nothing but gauge invariance. Actually, this should automatically lead 
in the quantum theory to the same renormahzation of the kinetic part and 
the interaction term in the action, Zkm = Zmt- It is not difficult to check 
that this is satisfied in our model, to the perturbative order in which we are 
working. From expressions (|27| ) and (^2]) we have 

Zk.n = l + |^(l-2.^){2F(l,l;l;^^)-7r^F(|||.2)}l + 0((4 

/ 9x f TT /I 1 1 „\ 8 , / 5 . 

Zi^t = 1 (i-2v^)\ F — ■,v^] - -v^F [2,2;-;v' 

167r2^ ^ I i; V2' 2' 2' / 3 V 2' 

which become identical after use again of the relations among hypergeometric 
functions. 
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The gauge symmetry 

A^ + d,,e{t,r) (37) 

is manifest at the classical level. Provided that it holds for the bare action 
(plD , it ensures that Zkm = Zmt- It is easily seen, by taking a gauge parameter 
9 independent of t, that it also enforces the same renormalization of the 
Fermi velocity v in the different terms of the action. As long as one is able 
to preserve the gauge invariance in the quantum theory, this turns out to 
guarantee the renormalizability to all orders of the perturbative expansion. 

The rest of the renormalization coefficients in (^) can be obtained from 
the following remark. The fact that Z^, has to equal Zi^t to all orders in 
perturbation theory implies, as in standard quantum electrodynamics, that 
Zg^y^ = 1. The renormalization coefficient of the electromagnetic field 
is determined from the polarization tensor, whose first perturbative 
contribution shown diagrammatically in figure 5 is a finite quantity 

n,.(fc) = ^4^?M.-^)|fc| /' dx ^x{l - X) (38) 

An V K JO 

We conclude, therefore, that 

ZT = Z;' = l + (e^) (39) 

remarking the absence of renormalization of the electric charge to first order 
in perturbation theory. 



4 Scaling of the electronic interaction (I). Non- 
relativistic regime. 

In this section we seek the connection of our model to the real physics of 
electrons in a two-dimensional layer. A single graphite sheet or any fullerene 
aggregate of sufficiently large size is the appropriate physical instance to look 
at. The first question concerns the magnitude of the parameters pertinent 
to the real layer. A standard value for the hopping parameter in graphite 
is t ~ 2.2eV . For the nearest neighbor separation in the lattice we can 
take a ^ 0.14 ■ 10~^m . From these values and restoring for a moment the 
dependence on c and h, the Fermi velocity in our model turns out to be 
V = 3ta/{2h) ~ 1.5 ■ lO^^c . This is certainly much smaller than the speed of 
light, for which to extract relevant quantities to the electron physics in the 
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layer we will have to take the limit f/c ^ and keep the leading order in 
previous formulas. 

For a clearer description of our approach to renormalization group we 
will quote the renormalization coefficients in this section as functions of the 
ultraviolet momentum cutoff A rather than in terms of the dimensional regu- 
larization poles At the one- loop level this just amounts to replace l/ehj 
log A in all formulas. The renormalization coefficient of the Fermi velocity 

becomes, for instance, in the nonrelativistic limit 

Zv = l-- logK + o[-] (40) 

It is not difficult to see that at each level of the perturbative expansion 
the leading order in the nonrelativistic approximation enters with a power 
(e^/f We define, therefore, a consistent perturbative expansion in the non- 
relativistic limit by taking t> — > (or, equivalently, c oo) and /v = const. 
This latter quantity becomes then the effective strength of the electronic in- 
teraction. In this limit we obtain from (EEi) 



Z^ = l + 0{^] (41) 




so that we do not find wavefunction renormalization at the one-loop level in 
the nonrelativistic limit. We recall also the former result f|39f) 



= Z;' = l + (e^) (42) 

We have a nonrelativistic quantum theory with no renormalization of the 
electric charge, but in which the effective strength of the interaction e^/ (Atcv) 
is renormalized. We rely on the original interpretation of renormalization 
adopted by Wilson to determine the scale dependence of the effective 
interaction. Since we have shown that the model is renormalizable, we may 
impose the independence of the renormalized theory on the momentum cutoff 
A. By means of renormalization group transformations we are able therefore 
to relate bare theories at different scales (different cutoff) which correspond 
to the same renormalized theory |T6|, p!7|. The essential information is encoded 



in the renormalization group equation, which for the electron propagator, for 
instance, states that under a change in the scale of the cutoff A 

(a A + p^^y^ e2)|- + (3,{v, e')^ - ^{v, e')) G{u;, k, A; v, e') = (43) 



14 



where 

Mv,e') = A^vn (44) 
Pe{v,e') = A^en (45) 
7(^,e2) = A-^logZ^ (46) 
The well-known solution to this equation takes the form 

G'(^,fc,pA;t;,e2) =exp|y''^7| G{uj,k,A;v,ff{p),elff{p)) (47) 

where the effective parameters are given by 

d 

d 



eff) 



P'o'p^^ff^P^ = -l3e{Veff,e^ff) (48) 

Increasing the scale p of the cutoff is equivalent to measure the observ- 
ables of the theory at large distance scale. From the one-loop order results 
of the previous section we see that e^jf is constant at this level, while v^ff 
is not. We expect, therefore, the effective coupling of the electronic interac- 
tion e^/(47rf) to have a nontrivial renormalization group flow in the infrared 
regime. Actually, from equations (HOf) and (|3) 



P^ =T^ + ^\^] ^9 

dp Vr IQTTVeff \KffJ 

SO that the asymptotic behaviour of the coupling is 

(A/Ao) = -\ + log— 50 

4vrfe// 47ruR mvR Aq / 

We have to bear in mind that, when integrating equation (l49|) , we are 
relying heavily on perturbation theory. If we recall our original estimate of 
the Fermi velocity f , we see that the coupling e^/ (47rf ) of the nonrelativistic 
theory is not a small parameter but, in fact, larger than one. Nevertheless, 
we know that when dealing with a renormalized quantum theory the state- 
ment that a coupling like e^/ (47rf ) is small or large applies to a particular 
energy scale. We then have to be very precise in establishing the regime in 
which the integration of the renormalization group equation becomes con- 
sistent. Let us remark that the solution (|50D of the renormalization group 
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equation is equivalent to performing the sum of leading logarithm terms of 
the perturbative expansion and curing, therefore, the problem of the infrared 
instability of the density of states addressed at the end of section 2. Taking 
the limit A — > oo in the bare theory corresponds to measure at large dis- 
tances (or small energies with respect to the Fermi level) and in this regime 
perturbation theory becomes more and more accurate. This is another in- 
stance of a physical system which has strong coupling at certain energies and 
an asymptotic regime at which perturbation theory becomes reliable, similar 
to what is believed to be true (as supported by experiment) in the quantum 
field theory of the strong interactions. We are providing the example of a 
nonrelativistic theory in which this also happens, exchanging the ultraviolet 
regime in the case of the strong interactions by the infrared regime here. 

As for the practical application of the flow (pO]), we remark that the 
nonrelativistic theory appears to be formally incomplete, since it does not 
give any hint about what the value of vr has to be. This is related to 
the unlimited growth of the Fermi velocity, and to the fact that, in the 
nonrelativistic approximation, there is no signal of a renormalization group 
fixed point. Later we will see that in the complete quantum field theory the 
fixed point is given by Vejf = c, but the flow in the graphite layer may not be 
necessarily controlled by this value of the renormalized coupling. We believe, 
anyhow, that the observation of the phenomenon of scaling in the asymptotic 
region should be possible in a single graphite sheet. 

Regarding the possible existence of a natural infrared cutoff, it can be 
shown that in our model the Coulomb interaction is not screened, in the 
usual metallic manner, by quantum effects. This is due to the vanishing of 
the density of states at the Fermi level. Thus, a finite screening length is 
only generated by giving to the system an excess of charge, as a consequence 
of doping. In order to study this effect, we may place the system away from 
half-filling by introducing a chemical potential /i for the conserved charge 

- = - lim Tr (T^(t,0W(t',0)) (51) 



The correct expression of the free electron propagator at finite n ispO| 



^ ^) .^0 = ^rzTTrf^ ~" ^\ '^'^ 

(52) 

Obviously, the /i-dependent part cannot change the ultraviolet behaviour 
of the loop-integrals and the renormalization coefficients of the quantum 
field theory are not modified by the nonvanishing chemical potential. This 
may lead, however, to finite corrections. The fermion loop diagram is now 
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nonvanishing, opposite to what happened before at /i = 0. The computation 
of (|5lD gives 



which, after multiplying by 2 (the two independent Fermi points) gives the 
correct number of particles. Regarding the polarization tensor, a nonvanish- 
ing value of noo(i^ = 0, — >• 0) signals the appearance of a photon effective 
mass. At small n one finds actually the estimate for the inverse of the screen- 
ing length 

2 

noo(o,fc^o)^^4 (54) 

which gives a measure of the screening effects away from half-filling. 

We now turn, finally, to the issue of whether the electron observables 
may get anomalous dimensions in the nonrelativistic limit. We remark that 
the absence of one-loop wavefunction renormalization in the limit v ^ is 
just an accident. The signal of wavefunction renormalization is given by an 
ultraviolet divergent self-energy contribution of the form 

S(cu,fc)^/(^^^7ou; + ... (55) 

To establish the degree of divergence of the local contributions one can safely 
implement the nonrelativistic limit in the photon propagator, which leads to 



the expression ([191) . The cj-dependence in the vertex diagram in figure 3 
can be absorbed then by a simple change of variables in the loop integral, 
so that the divergent part of the diagram cannot bear any dependence on 
it. At the two- loop level, the same independence of the local contributions 
on u holds except for the vertex diagrams in figure 6. In these diagrams 
the loop integrals overlap in such a way that the previous argument does 
not apply. By using the nonrelativistic approximation (p!9|) and after a very 



lengthy calculation, one obtains from the sum of the two-loop diagrams 

^* = l-^^/o^7A + Oie') (56) 



where ci is a nonvanishing quantity, ci ~ 5.49 ■ 10~^. We see therefore 
that the electron observables acquire anomalous dimensions in the infrared 
regime. According to equation (|T7|), the electron propagator transforms un- 
der a change of scale p by the factor 



eff 
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We have therefore an electronic system which deviates shghtly from Fermi 
hquid behaviour. In particular, we have seen that the vanishing of the density 
of states at the Fermi level suppresses the usual screening of the Coulomb 
interaction. This property is maintained in the quantum theory, at least at 
the one-loop level, since the Fermi level is not shifted by quantum correc- 
tions. Also a gap does not open at the Fermi points to the one-loop order (as 
it would be if an effective mass for the spinor had been generated in pertur- 
bation theory). The system is not an insulator, though it does not exhibit 
usual metallic properties either. 

5 Scaling of the electronic interaction (II). 
Non- Fermi liquid fixed point. 

In the previous section we have mentioned that the nonrelativistic theory 
is formally incomplete since it does not know of any bound in the effective 
Fermi velocity v^ff of electrons. As this parameter grows steadily at large 
distances the nonrelativistic approximation becomes less accurate, up to a 
point in which the power series expansion in v/c looses its meaning. In this 
relativistic regime, we have to use the full expression of the beta function 
Pv{v,e^) to analyze the behaviour of Vgff. In particular, we are interested 
in knowing if a possible fixed point of the renormalization group flow exists, 
characterized by the condition 

P.{v,e^)=0 (58) 

By quoting (pOf ) in terms of elementary functions, this equation reads, at the 
one-loop level, 

1 1 - 2i;2 + 4:v^ 1 - Av^ ^ 

n arccos v + ^ = 59) 

For positive v, the only solution to (^) is f = 1 [], which in our units means 
V equal to c. Thus, the full quantum field theory has a built-in upper bound 
of the renormalization group flow, and the speed of light turns out to be the 
fixed point for the Fermi velocity of electrons. 

In fact, it is not difficult to check that the one- loop beta function has a 
simple zero at f = 1, and to compute the fiow in the neighborhood of the 
fixed point. The renormalization group equation for Vefj becomes 

op Vr 71^ 



''This is, obviously, a statement regarding the hmit value of the left hand side of 
as w ^ 1. 
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The asymptotic behaviour of the solution is, taking for vr the fixed point 
value, 



In figure 7 we show the plot of this scaling behaviour, which is practically 
indistinguishable in the range considered from that obtained with the non- 
relativistic approximation. 

Recalling the vanishing of the electric charge beta function /3e(f,e^) at 
the one-loop level, we may assert that to this perturbative order we have 
identified an infrared fixed point given by Vejf = c and arbitrary (though 
small) value of Ce//. The one- loop result is not satisfactory, however, since 
the independence of Ce// on the flow may be destabilized by higher- loop 
effects. In the renormalization group approach to Fermi liquid theory the 
consideration of the one-loop graph in figure 5 is enough since it encodes, via 
the renormalization group, information of the series of ladder diagrams. In 
the usual analysis of the many-body theory with an extended Fermi surface 
these are which dominate in the infrared regime. Because of the existence 
of Fermi points in the honeycomb lattice at half-filling, though, there is no 
argument here to conclude the relevance of one set of diagrams against the 
rest. In our quantum field theory description it seems that no diagram can be 
discarded a priori. It is also for this reason that there are diagrams which are 
irrelevant in Fermi liquid theory but give rise in our model to wavefunction 
renormalization and anomalous dimensions in the infrared regime. 

Regarding the renormalization of the electric charge, we have analyzed 
the polarization tensor n^,^ to second order in perturbation theory. As stated 
above, the renormalization of the photon field Afj_ is directly related by gauge 
invariance to the renormalization of the charge e. Computing in the frame- 
work of dimensional regularization we have not found any divergent contri- 
bution of order to 11^,^. The reason is the same operating in the one- loop 
polarization tensor. A naive dimensional analysis of the diagram shows that 
it should diverge linearly in momentum space. In the dimensional regulariza- 
tion scheme this leads to the appearance of gamma functions at half- integer 
values of the argument, instead of poles at d ^ 3. We beheve that this 
scheme is the best suited for the regularization of the gauge theory — it cer- 
tainly is not going to break the gauge invariance in any event. If, however, 
one insists on understanding the regularization in terms of a cutoff depen- 
dence, it can be checked that this does not lead either to a renormalization 
of the scale of the field A^. One would have to look for divergent contribu- 
tions of log A type accompanying the gauge invariant expression (^) for 
Hf^i^. These are not present at the two-loop order and one may conjecture, 
following the above dimensional argument, that they do not appear to any 




(61) 
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order in perturbation theory. 

At or near the fixed point, the two-dimensional electron liquid we are 
describing is definitely not a Fermi liquid. The main reason for this assertion 
is the anomalous dimension acquired by the electron field which may be 
obtained from (27) evaluated now at f = 1. At the fixed point the behaviour 
of the Green function under a change of scale is 

G{uj,k)\^ = A''<l>{uj,k) (62) 

or, equivalently, after introducing the precise value of 7 

G{puj, pk) = p-i+'^'/(^'"')G'(cu, k) (63) 

We arrive in this way to solve the problem of the proliferation of ^05^0; terms 
in the perturbative computation of the density of states. As is well-known, 
the renormalization group approach succeeds in summing up all these increas- 
ingly divergent contributions, in the form of an anomalous dependence on 
the electron frequency. We have found that this sum is consistently achieved 
in the relativistic regime, where the fixed point lies and the electron field has 
a simple, anomalous scaling behaviour. 



6 Summary and conclusions 

We have studied a system of electrons in the two-dimensional honeycomb 
lattice interacting via a Coulomb potential. The fact that the Fermi sur- 
face of this lattice at half-filling consists of two isolated points renders the 
problem similar to the d = 1 case and allows the formulation of the low 
energy effective theory as a standard quantum field theory. We have been 
able to investigate the infrared behavior of the system by using renormaliza- 
tion group techniques and have seen that, indeed, it resembles more that of 
one-dimensional electrons than that of a two-dimensional Fermi liquid. 

Before entering into the analysis of the phenomenological consequences 
of our work, there are a few points that we would like to emphasize. 

Let us first comment on the use made through the paper of the electro- 
magnetic field propagator. This use is unconventional both from the con- 
densed matter point of view where the photon is mostly absent, and from 
the quantum field theory since, as we will strength now, we are not dealing 
with standard massless QED. 

Expression (0) shows the usual three-dimensional photon of QED. This 
is what we use in our relativistic computations. The constraint that, in the 
electronic interactions, the two points lay on the lattice plane, allows us to 
integrate out the perpendicular dimension. The resulting expression read 
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as a propagator on the plane certainly differs from the standard propagator 
in QED3. In particular its infrared behavior is quite different, it goes like 
1/r instead of the standard log r. This is the origin of, for instance, the 
finiteness of the photon self-energy up to the two-loop level and, hence, of 
the nonrenormalization of the electric charge. 

Let us now comment on the comparison between the instantaneous versus 
the retarded propagator. The effect of retardation on the infrared behavior of 
nonlocal quantities has been exemplified with the computation of the quan- 
tum corrections to the density of states at the end of section 2. In the theory 
given by the Coulomb interaction (|l^), no one- loop corrections were found 
to the free value of the density of states. In the theory with the retarded 
interaction our RG analysis shows that, despite the infrared divergence found 
at the one-loop level, the approach is consistent in the sense that the pertur- 
bative expansion of the nonrelativistic approximation does not have to make 
sense in the infrared when the Fermi velocity is still far from the fixed-point 
value. Hence, we must resort to the full relativistic theory when computing 
the infrared behavior of nonlocal observables. There, the RG allows to inter- 
pret the one-loop infrared instability shown in section 2 as the manifestation, 
in perturbation theory, of the anomalous dimension. 

We have then found that our original system flows in the infrared towards 
a fixed point of the renormalization group that is a relativistic quantum field 
theory with anomalous dimension of the fermionic propagator. 

Our analysis has important consequences in condensed matter physics. 
One of the more interesting is the anomalous screening of the Coulomb in- 
teractions of the model. In a Fermi liquid, formation of virtual electron-hole 
pairs (i. e., quantum corrections to the photon self-energy) generates a new 
scale, the screening length. This is the basis of the absence of long-range 
interactions in the ultralocal approach to interacting electrons (Hubbard 
model). What we have seen is that, due to the vanishing of the density 
of states at the Fermi level and to the absence of renormalization of the 



photon self-energy, there is no metallic screening in our materials [0. The 
robustness of this result stands on the non-renormalization of the chemical 
potential (the Fermi level remains a set of two points after the corrections 
due to the interactions). Moreover we have seen that, although the low en- 
ergy density of states goes to zero, no finite gap appears in the spectrum of 
charge excitations. Our results rule out the validity of descriptions in terms 
of on-site interactions only. Screening effects are not strong enough to reduce 
the couplings to simple on-site terms. 

The one-particle properties are very different from those of a conventional 
Fermi system. The wavefunction renormalization leads to the disappearance 
of the quasiparticle peak at low energies. Thus, no effective description in 
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terms of independent, and coherent, quasiparticles is possible. The situation 
resembles the ID Luttinger liquid, where any amount of electron-electron 
interaction suffices to wipe out the quasiparticle pole, leading to incoherent 
one-particle excitations. 

Finally, the existence of nontrivial scaling laws implies that anomalous 
exponents should show up in a variety of properties, like the specific heat, 
susceptibility, etc. The possibility of having an experimental test of these 
results remains open. 

To conclude this summary and to point out possible future work, we 
would like to remark the interest of investigating the connections of our work 
with the infrared catastrophe mentioned in the introduction as well as with 



other unconventional Fermi liquid behavior discussed in the literature 



On this respect, this work can be seen as a rigorous realization of the ideas 
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Figure Captions 



Figure 1: The planar honeycomb lattice. 

Figure 2: Representation in {E, k) space of the lower branch of the 
electronic dispersion relation (in units t = a = 1). The cusps appear at the 
six corners of the first Brillouin zone. 

Figure 3: One-loop contribution to the electron self-energy. 

Figure 4: First quantum corrections to the interaction vertex. 

Figure 5: One-loop polarization tensor. 

Figure 6: Two-loop diagrams contributing to wavefunction renormaliza- 
tion. 

Figure 7: Plot of the effective Fermi velocity Ve/f (in units c—1) versus 
length scale L/Lq. 
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